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The spin-polarized transport in a single-molecule magnet Fe4 sandwiched between two gold electrodes is 
studied, using nonequilibrium Green's functions in combination with the density-functional theory. We 
predict that the device possesses spin filter effect (SFE), spin valve effect (SVE), and negative differential 
resistance (NDR) behavior. Moreover, we also find that the appropriate chemical ligand, coupling the single 
molecule to leads, is a key factor for manipulating spin-dependent transport. The device containing the 
methyl ligand behaves as a nearly perfect spin filter with efficiency approaching 100%, and the transport is 
dominated by transmission through the Fe4 metal center. However, in the case of phenyl ligand, the spin 
filter effect seems to be reduced, but the spin valve effect is significantly enhanced with a large 
magnetoresistance ratio, reaching 1800%. This maybe attributed to the blocking effect of the phenyl ligands 
in mediating transport. Our findings suggest that such a multifunctional molecular device, possessing SVE, 
NDR and high SEE simultaneously, would be an excellent candidate for spintronics of molecular devices. 

Molecular spintronics using molecules as spin transport elements has attracted intensive interest due to 
obtained various functionality, and molecular devices such as spin valve, single-molecule transistors, 
switch and diode, etc^ l Conventionally, the magnetic molecule junction (MMJ) can be achieved by 
placing a nonmagnetic molecular bridge between two ferromagnetic electrodes^'^. In the rather newer field, MMJ 
can also be achieved by coupling a magnetic molecule to two nonmagnetic probes^"^°. In this approach, single- 
molecule magnets (SMMs) are ideal building blocks to construct MMJ because of their independent magnetically 
functional units^^"^^, where the inner magnetic core in SMMs is surrounded by organic ligands and the interaction 
between magnetic cores of neighbor molecule is very weak, leading to a large- spin ground state. Unlike above 
mentioned conventional MMJ depending on reasonable magnetic configuration of two ferromagnetic leads^'^^, 
the origin of magnetic behavior in this MMJ consisting of a single-molecule magnet and two nonmagnetic 
electrodes is the intrinsic spin of the center magnetic molecule, therefore avoiding the great difficulty of the spin 
injection from ferromagnetic electrodes into a nonmagnetic center in real applications. And so the unique feature 
of MMJ containing SMM would offer much greater potential for high-density information storage devices^^. 

Although there have been some successful efforts to organize SMMs on various substrates, the observation of 
magnetic hysteresis on individual molecules organized on surfaces^^'^^, that is a necessary step to develop molecu- 
lar memory arrays, had not been reported for monolayers of SMMs when wired to metallic surfaces until Matteo 
Mannini and co-workers found that tailor-made tetranuclear Fe4 complexes retain their magnetic properties at 
gold surfaces This encouraging result naturally excites our attention to investigate quantum transport through 
this high-performance magnetic molecule Fe4. Several recent works mainly focused on the charge transport and 
the magnetic anisotropy in a single- electron transistor containing this magnetic molecule^^"^". However, to the 
best of our knowledge, the spin-polarized transport on this magnetic molecule has not been reported both 
theoretically and experimentally so far. Here, we carry out first-principles transport calculations through the 
single magnetic molecule Fe4 linked to two gold nanowires. Furthermore, we firstly explore the potential of a spin- 
polarized current through this magnetic molecule to read and write spin-based information. We predict that this 
device has a nearly perfect spin filter effect with efficiency approaching 100%. Moreover, it also shows spin valve 
behavior and negative differential resistance effect. In addition, we also find the ligand group coupling the 
molecule to leads is a key factor for manipulating spin dependent transport. These findings strongly suggest that 
such a multifunctional molecular device would be candidate for molecular spintronics. 
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Results 

The structure of Fe4 molecule is shown in Fig. 1(a), where the metallic 
cores are composed of a triangle of iron ions and an inner ion^\ Two 
derivatives are studied: Fe4Me functionalized with tripodal ligands 
methyl-C(CH20H)3 (See Fig. 1(b)), and Fe4Ph functionalized with 
phenyl-C(CH20H)3 (See Fig. 1(c)). The two ligands have been 
shown to be advantageous for deposition of Fe4 molecule on gold 
surfaces^^'^^. The device we considered consists of two semi-infmite 
gold electrodes sandwiching each Fe4 derivative via thiol end groups, 
as shown in Fig. 1(b) and (c). The gold electrodes extend to z = ± 
along the transport direction, and extend with large vacuum layer in 
transverse x and y directions in order to suppress the interaction 
between transport system and its mirror images. 

First, we study the transport properties of Fe4Me derivative, where 
the Fe4 core is connected to the gold nanoelectrodes via methyl 
ligands terminated with thiol groups as shown in Fig. 1(b). Two 
magnetic structures are considered for the spin- dependent transport 
properties. One is antiparallel configuration (APC) as depicted in the 
left panel of Fig. 1(a), where the outer individual spins interact anti- 
ferromagnetically with the inner one, to give a S = 5 ground spin 
state. Another is parallel configuration (PC) as seen in the right panel 
of Fig. 1(a), where all the four iron cores have the parallel spin 
direction, to give a S = 10 high-spin state. A magnetic field can switch 
the Fe4 molecule from one state to the other, and the further calcula- 
tions indicates that the magnitude of magnetic field needed to switch 
the molecule from the antiferromagnetic state to the ferromagnetic 
state is about 32 T, which is comparable with the results obtained by 
Baranger et af. Fig. 2(a) and 2(b) show the current-voltage (I-V) 
curves of the magnetic molecular junction for APC and PC, respect- 
ively. The total current for APC and PC are rather different, resulting 
in a large magnetoresistance (MR) ratio defined as MR = (Iapc ~ 
Ipc)/Ipc as shown in Fig. 2(c), where Iapc and Ipc are the total 
currents for APC and PC related to the same applied voltage, respect- 
ively. At zero bias when all currents vanish, MR ~ 75% is obtained by 
using equilibrium conductance. This behavior shows that the device 



(a) APC PC 




Figure 1 | The schematic illustration of the SMM device, (a) A view of the 
magnetic core structure of Fe4 with antiparallel configuration (APC) and 
parallel configuration (PC), (b) Fe4Me functionalized with methyl ligands, 
and (c) Fe4Ph functionalized with phenyl ligand are respectively bonded to 
the two semi-infinite Au leads via thiol end groups. H atom is not shown 
for the sake of clarity. 



can also be used as an organic spin valve. The calculated zero-bias 
conductance for APC remains essentially the same magnitude as the 
experimental results Furthermore, Fig. 2(b) displays that the total 
currents for PC are only from the spin-down channel, the spin-up 
current is completely suppressed, leading to a good spin filter effect. 
Also, a spin-filter diagram (see Fig. 1 in Ref. 22) can be used to 
describe how the spin current in Fig. 2 (b) tunnels through the device. 
The spin filter efficiency defined as = (J| — I|)/(Ij + I|) for PC is 
higher than 99% in the whole bias region as shown in Fig. 2(d), 
indicating a perfect SFE. In addition, as seen from Fig. 2(a) and 
2(b), the total currents for APC and PC increase nonlinearly with 
increasing bias voltage, and reach a maximum around 5.6 mV. 
Afterward, the total currents decline when the applied bias voltage 
is increased to larger values, leading to a negative differential resist- 
ance (NDR) behavior. 

It is well known that the ligand coupling the molecule to the leads 
is a crucial issue for the preparation of single-molecule electronic 
device. Here, we also investigate the transport properties of Fe4Ph 
derivative by substituting the methyl groups in Fe4Me with phenyl 
ligands (See Fig. 1(c)). The calculated results for spin parallel con- 
figuration and antiparallel configuration are given in Fig. 3. 
Compared with the above-mentioned Fe4Me in Fig. 2, the transport 
property of Fe4Ph derivative shows two main differences: (1) The 
strong spin polarization appears in APC rather than in PC, different 
from the situation in Fe4Me depicted in Fig. 2(a) and 2(b). 
Furthermore, we discover that the spin filter efficiency for PC in 
Fe4Me (See Fig. 2(d)), is higher than that for APC in Fe4Ph (See 
Fig. 3(d)), which indicates the methyl groups bridging the Fe4 mole- 
cule and electrodes are more advantageous for the spin filter effect 
than the phenyl ligands. (2) In PC of Fe4Ph system, both the spin-up 
and the spin- down currents are limited, and thus the total current for 
APC is much larger than that for PC, exhibiting a substantial mag- 
netoresistance with a maximum 1800% (See Fig. 3(c)), which is 
significantly larger than that observed in Fe4Me. This suggests that 
the phenyl ligand is highly advantageous for the spin valve effect. 

Discussion 

To further explore the underlying mechanisms of MR and SFE in 
Fe4Me device. Fig. 4(a) and 4(b) present the spin-resolved transmis- 
sion coefficient T(£) and the spatial distribution of the molecular 
projected self- consistent Hamiltonian states (MPSHs) of the HOMO 
(highest occupied molecular orbital) and LUMO (lowest unoccupied 
molecular orbital) at zero bias. For APC, the T(£) of Fig. 4(a) shows 
that the total current comes from both the spin-up and the spin- 
down transmission peaks around the Fermi level. Furthermore, 
MPSHs in Fig. 4(a) reveal that the outer three iron ions provide 
the spin-down transmission channels, and the centered iron ion 
contributes to the spin-up channel. Whereas, for PC in Fig. 4(b), 
T(£) and MPSHs indicate that the resonant transmission peaks are 
only from the spin-down electrons of the four iron ions, the spin-up 
T(£) is essentially zero around the Fermi energy, which suggests that 
only spins parallel to the molecular magnetization can flow through 
the SMM, exhibiting a very high spin polarization with efficiency 
approaching 100%. Simultaneously, this also results in that the total 
current in APC is significantly larger than that in PC. And so the 
device can be used for a spin valve by switching between APC and PC 
magnetization configurations under an applied magnetic field. 

In order to explore the physical origin of the NDR phenomenon, 
we investigate the spin-resolved transmission properties around the 
Fermi level. Since the NDR behaviors are found for APC and PC 
almost at the same bias region, we only give the transmission coef- 
ficient of PC system at 6.8mV and 5.6mV, as shown in Fig. 5(a) and 
5(b), respectively. One can see that the main transmission peak, 
corresponding to the conduction orbital LUMO (lowest unoccupied 
molecular orbital) above the Fermi energy, is responsible for the 
NDR behavior. The LUMO peak originating mainly from the 
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Figure 2 | The transport properties of Fe4Me. (a) and (b) Current-voltage characteristics for magnetic antiparallel configuration and parallel 
configuration, respectively.Total means the sum of the spin-up and spin-down currents for each magnetic configuration, (c) Magnetoresistance ratio, (d) 
Spin filter efficiency for parallel configuration. 



orbitals of four Fe cores is higher in Fig. 5(b) than in Fig. 5(a). In deciding which parts of the chemical structure should be responsible 

addition, we notice that a short transmission peak appears at 0.08 eV for the transmission. First, we study the local contributions to trans - 

in Fig. 5(b), which arises from the hybrid state between the molecule mission for the Fe4Me system. For APC, This is, in fact, exactly what 

and the electrode surface, whereas, the peak is not found in Fig. 5(a). is seen in Fig. 6(a) and 6(b), the spin down density is localized on the 

Therefore, the hybrid state is another reason for the NDR behavior. outer three iron ions, and likewise, the spin up electron density is 

Now, let us turn attention to the influence of the ligand on the localized on the centered iron ion with only a little distribution on the 

transport properties. The local density of states (LDOS), as depicted outer three iron ions. The fact indicates that the total transmission in 

in Fig. 6, provides valuable information, which can be helpful in APC is from the sum of contributions of the spin-up and spin-down 




0246802468 
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Figure 3 | The transport properties of Fe4Ph. (a) and (b) Current-voltage characteristics for magnetic antiparallel configuration and parallel 
configuration, respectively. (c) Magnetoresistance ratio, (d) Spin filter efficiency for parallel configuration. 
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Figure 4 | Spin-polarized transmission coefficient and MPSH states of 
Fe4Me at zero bias (a) for antiparallel configuration and (b) for parallel 
configuration. The average Fermi level is set to zero. 

electrons. For PC, Fig. 6(c) and 6(d) illustrate that only one spin 
direction (i.e. spin up) of the conducting electrons is transmitted, 
leading to a nearly perfect spin filter effect, simultaneously, explain- 
ing the occurrence of spin-valve behavior. These behaviors are in line 
with the results from the energy spectrum of the molecular projected 
self-consistent Hamiltonian (MPSH) at zero bias (Supporting Infor- 
mation, Figure SI). In general, the transport in Fe4Me system 
whether for APC or PC is dominated by transmission through the 
Fe4 metal center, less sensitive to the methyl ligand. However, the case 
is different in Fe4Ph device, except for Fig. 6(e), where the metal centers 
are still the main contributor for the transmission. From Fig. 6(f) to 
6(h), and in particular for Fig. 6(f), the through-metal contribution is 
well below the through-ligand contributions. Furthermore, for all 
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Figure 5 | The bias-dependent transmission coefficient of Fe4Me in 
parallel configuration (a) for 6.8mV and (b) for 5.6m V. The spin-up 
transmission coefficient is almost zero. 
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Figure 6 | Surfaces of constant LDOS at the Fermi energy, (a/b) and (c/d) 
for device Fe4Me depict density of states in antiparallel configuration and 
parallel configuration, respectively. (e/f) and (g/h) are similar to (a/b) and 
(c/d), but for Fe4Ph. The inset on the right side of (f) shows the sites of 
carbon atoms in phenyl ring. Green (blue) represents spin down (spin up) 
density of states. 

three cases, it is found that the densities are localized on 1, 3, 5 sites 
in the left phenyl ring, forming the conjugated n bonding orbital, and 
likewise, the densities are localized on 2, 4, 6 sites in the right phenyl 
ring, forming the tt* antibonding orbital. Since the antibonding 
orbital has higher energy than the bonding orbital, it is difficult that 
either spin from the left phenyl ring penetrates into the right one, and 
hence, the conductance is low. This implies that the blocking effect of 
the phenyl ligands on the electronic transport in Fig. 6(f) to 6(h), is 
mainly responsible for the large spin valve and spin filer effects in 
Fe4Ph system. 

Methods 

The spin-dependent electron transport calculations are performed using a recently 
developed package Atomistix Toolkit^^"^^, which combines the nonequilibrium 
Green's-function (NEGF) with a self- consistent density- functional theory (DFT)^*^'^^. 
The generalized gradient approximation (GGA) in Perdew-Burke-Ernzerhof (PBE) 
format^*^ is adopted for the exchange correlation potential in DFT. Further details of 
the method were presented in Ref 29. The wave functions are expanded on a 
numerical basis set of double zeta-polarized basis. The mesh cut-off energy is 150 Ry 
and the convergent results are achieved by using the k-points sampling 100 in 
transport z direction. This procedure is iterated until it reaches self- consistency when 
the convergence criterion 10"^ for the total energy is achieved to present an accurate 
charge density. The geometry of the scattering region is relaxed when the force 
tolerance of 0.05 eV/A is achieved, where the outmost layers of Au atoms are fixed at 
their bulk positions. After the self- consistent NEGF-DFT procedure is converged, we 
obtain spin-polarized current using the Landauer-Buttiker formula: 



la-- 



j-a 



TAB, n) IfdE - III) -fR(E - fiR)] dE, 



where fi R is the Fermi-Dirac distribution and fii^R are the chemical potentials of the 
left electrode (L) and the right electrode (R), respectively. T(j{E,Vb) is the spin- 
dependent transmission coefficient at the energy E with bias voltage Vb- 

Finally, we attempt to check the DFT results by using a different exchange cor- 
relation functional. For comparison, we also apply the local density approximation 
(LDA) in Perdew-Zunger parametrization as the exchange correlation functional, 
and investigate the transport properties of Fe4Me device as an example. The calcu- 
lated I- V curves with the LDA functional are provided in the Supporting Information 
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(Figure S2), where one can observe the NDR peak whether for APC or PC are shifted 
to higher energies (compare the NDR peaks with the GGA functional in Fig. 2(a) and 
Fig. 2(b)). Simultaneously, one can see that the current value (a maximum of 18nA) 
with the GGA functional is larger than that (a maximum of 12nA) with the LDA. 
Accordingly, a MR maximum of —80% with GGA is reduced to —60% with LDA. 
However, the overall shape of these transmission curves using LDA remains essen- 
tially the same as those using GGA. Furthermore, a MR maximum of —85% obtained 
with LDA is also in line with results from GGA, though the calculated current with 
LDA is larger than that with GGA. This indicates that changing the functional will 
result in quantitative changes in the transmission, but not affect our qualitative 
conclusions, such as the nature of the NDR, spin valve, and spin filter effects. 
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